Spinor BECs in a double-well: population transfer and Josephson oscillations 
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The dynamics of an F = 1 spinor condensate in a two- well potential is studied within the frame- 
work of the Gross-Pitaevskii equation. We derive two-mode equations relating the population im- 
balances, the phase differences among the condensates at each side of the barrier and the time 
evolution of the different Zeeman populations for the case of small population imbalances. The 
case of zero total magnetization is scrutinized in this limit demonstrating the ability of a two mode 
analysis to describe to a large extent the dynamics observed in the Gross-Pitaevskii equations. It is 
also demonstrated that the time evolution of the different total populations fully decouples from the 
Josephson tunneling phenomena. All the relevant time scales are clearly identified with microscopic 
properties of the atom-atom interactions. 
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The recent experimental work on optically trapped 
spinor condensates has broadened the frontiers of con- 
fined Bose Einstein condensates (BECs) P, 0. There, 
a transfer of population between the different Zeeman 
components of a spinor BEC was observed providing a 
clear signal of the spin-dependent interatomic interac- 
tions. This experiment quickly connected the field of 
cold-atoms to a large variety of problems in quantum 
magnetism, mostly related to magnetic ordering and spin 
dynamics [3] . From a different point of view, it can be re- 
garded as the first case of ternary mixture of BECs with 
population exchange among the three components, a nice 
example of coupled multicomponent quantum gases. 

At the same time, the fast tunneling of atoms through 
potential barriers driven by imbalanced populations at 
each side of the barrier was shown experimentally, short 
after it was observed in optical lattices [3], providing 
the first confirmation of Josephson tunneling of atoms 
in BECs [5:]. From the BEC point of view, Josephson 
tunneling through the potential barrier produces a weak 
coupling between the BECs at each side of the trap, pre- 
senting a coupled-oscillator behavior in the appropriate 
variables Q. 

Both effects have been considered by several theoretical 
groups in quite different contexts, e.g. d, 0, H, H, [13, [Hlj 
In this letter we address both aspects together: the 
Josephson oscillations and the transfer of populations. 
Providing a direct connection between the different time 
scales and the microscopic properties of the interatomic 
interactions. 

To fix the conditions we consider the simplest scenario, 
which already contains relevant physics. We restrict our 
analysis to the case of small population imbalance of all 
the Zeeman sublevels and small initial phase difference 
between the same component at both sides of the trap. 
Our main tool are the mean- field Gross-Pitaevskii (GP) 
equations for spinor BECs On top of numerical 

solutions to the GP equations we illustrate the physics 
emerging in this fairly complex situation by deriving a 
two-mode description of the problem. Within this two 



mode description it is easy to show that for small pop- 
ulation imbalances and phase differences the population 
transfer dynamics fully decouples from the tunneling phe- 
nomena. 

In the mean field approximation the dynamics of the 
vector order parameter 5* = (^'_i, \l/o, vE'i) representing 
the F — 1 spinor condensate is given by jlOj . 

ihdt-^±i = [T^s + C2(n±i +no -n^i)]*±i + C2*g*;,i 
ihdt-^a = [7^s + C2(ni+n_i)]«'o + C22'J'i«'S'J'-i,(l) 

where. Us = - + V + Cgn, n„(f, t) = |^'„,(f, i)|2, 
n{r,t) = X^m "mC^' ^'^'i = 0,±1. The population 
of each hyperfine sublevel is Nm (t) = J df rim {r, t) . Due 
to the last term in the r.h.s of Eqs. ([T]), the population of 
each Zeeman sublevel is not conserved. The couplings are 
Co = AiTti^{aQ + 2a2)l{iM) and C2 = 47rft2(a2-ao)/(3M), 
where oq and a2 are the scattering lengths describing 
binary elastic collisions in the channels of total spin 
and 2, respectively. Their values for *^Rb are ap = 
101. Sas and 02 = 100.4aB[il, which yield C2 < 0, 
thus producing a ferromagnetic-like behavior. The to- 
tal number of atoms in the system and total magneti- 
zation are conserved quantities, N ~ J df n{f,t) and 
M = J df[n+i{f,t) - n_i(r,t)] . 

We consider a setup similar to that described in Ref. Q 
but with two important differences: the total number of 
atoms and the barrier height. In our case the number of 
atoms is larger, N — 15000, in order to enhance popula- 
tion transfer effects. We use the same kind of double-well 
potential but with a higher barrier and a tighter con- 
finement in the x direction to ensure a clear Josephson 
tunneling situation. The potential then reads, 

Vir) = Y(t^'x2+w2y2+^2^2)^^^^Qg2(^^/^^) 

with tjj, = 27r X 100 Hz, ujy = 27r x 66 Hz, iv^ — 2tt x 90 
Hz, qo = 5.2/im, Vq = 3500 h Hz and M is the mass of the 
atoms. As in the experiment [5] we assume that the dy- 
namics takes place essentially on the x axis. Then, defin- 
ing uj±_ — ^ujzUjy the coupling constants can be rescaled 
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FIG. 1: (a) (above) 4>g and "3?ist. (below) Potential in the x 
direction together with one of the initial population profiles 
used in the simulations (arbitrary units). The solid (black) 
curves show the evolution of the total population (b) and of 
the population imbalances and phase differences (c) corre- 
sponding to simulation I of Tab. |T] The dashed (red) lines 
depict the two-mode calculation. 



by a factor l/(27ra^), with a± the transverse oscillator 
length [l4f|. and the dynamical equations transform to 
one-dimensional ones for a symmetric double-well. 

The numerical simulations of Eqs. (jT]) are performed 
in the following way. First using an imaginary time evo- 
lution method we compute the ground, ^gs = and 
first excited state, ^ist = $- of a scalar BEC, C2 — 0, 
under the same conditions. Then, given initial popula- 
tion imbalances for all the components, wc build t = 
wave functions by the appropriate linear combinations of 
and We study the time evolution of the sys- 

tem by means of the split operator method. In Fig. [T] 
we depict $gSj ^ist and the potential in the x direction 
together with one of the initial density profiles used in 
the simulations. 

To characterize the Josephson dynamics we de- 
fine for each component the population imbalance, 
Zm = (Nm^L - Nm,R)/Nm, and the phase differ- 
ence, l5(/)„ = (f)m,R ~ (l)m,L- WhcrC, Nm,L{t) = 
/-ooda^ /"co/-°ldydzn„(f,i), Nm,R{t) = ^mW " 
Nm,L{t) and 4>m,R(L) are the space average of the phase 
of 5'm(r, t) at each side of the barrier. Lets us emphasize 
that the phase of {r, t) is almost spatially constant at 
each side of the trap during the GP simulations. 

To better understand the dynamical content of the 
GP results, Eqs. H]), we derive a two mode approx- 
imation for the GP description of spinor BEC based 
on the following ansatz, ^'^(r, t) = '^m,L{t)^m,L{r) + 
*m,_R(t)^'m,fl(^) , where the modes $m,L(fl) are mostly 
localized at the left (right) side of the trap. These modes 
can be built from <f>± obtained from the GP equations, 

*m,L = + <^ni.R = -^{^^n,+ - $rn.-) 

with ^m,±{T) = ±<I>m^±(— r). The complex compo- 
nents are normalized as, ^m.L = \/Nm,L{t)e^'^'"-^ and 
'i'm,R = \/Nm^R{t)e"^"'-'^. Therefore, for this case, 

Zm = {Nm,L - Nm,R)/Nm, and S(f>m = (t^m.R ~ <p7n,L Vm. 

As a first step we consider the so-called standard two- 
mode, which implies that all the overlapping integrals 
involving products of L and R modes of any two compo- 
nents are neglected. This approximation is expected to 



yield essentially the correct physics as it was for the case 
of the scalar condensate or binary mixtures 0, • 

We take the following assumptions: zero total magne- 
tization, small imbalances, Zm, small phase differences, 
6(l)m, and small 6(j)L{R) = 20o,L(fl) - 4>+i,l{r) - 4>~iMR)- 
Also we take equal modes for the three states {'^l{r) = 
^m.L{R) Vm), which is fully justified at zero magnetic 
field (it corresponds to the single mode approximation at 
each well). 

In such conditions one can prove that the total popula- 
tion of the different components, Nm(t), fully decouples 
from the Josephson tunneling dynamics. The time evo- 
lution of A^o is given by. 



No{t) 



-WiNo{t){N - NoimNoit) - A/2) (2) 



with hU2 — C2 J df4'L(R)i^' where one can indistinctly 
use the left or the right mode. The other two follow: 
N±i{t) = (A - No{t))/2. If No{t) ~ A'/2, the behav- 
ior of A^o becomes sinusoidal, No{t) = A/2 + {No{0) - 
A/2) cos(tt'Ti), where we have defined the "population 
transfer frequency", ojt = NU2- Eq. ([2]) gives an excel- 
lent agreement compared to the full GP results. With 
the considered conditions no damping is observed in the 
evolution of the populations. 

The system of equations governing the dynamics of the 
population imbalances, Zm, and phase differences, 6(j)m, 
reads: 

Z±l = -LUrS(t)±l - (Ao/2) U2i5<j) + Z±lA(^) , 
Zq = —UJrS(f>0 + ^ U2{S(j) + ZoAcj)) , 

4±i = U{Nz±i + NoZo) + U'Nz^i 
No 

+UJrZ±i + U2 — {2zo - Z±i + Zzpi) , 

S<i>o = {U + U2)N{z^i+z+i) + UoNoZo + uJrZQ 
= 8{No-N/2)U2, (3) 

where Scfi — S(f>L — S(f)R, Acj) = 54'l + 54>r, N = N+i = 
N^i = {N- No)/2, hUo = Co / drcbl^j^^ir), U = Uo + 

U2, U' ^Uo- U2, K ^ - f dr{hy{2M)V<^L ■ + 
V ^r), and ojr = 2K/h, is the Rabi frequency. 

From the ground and first excited states of the sys- 
tem computed numerically, see Fig. [U we build the left 
and right modes as explained above and compute the mi- 
croscopic parameters entering in the two-mode descrip- 
tion. The resulting values are: ujr = 0.00386 KHz, 
NUq = 26.604 KHz and NU2 = 0.12366 KHz. This 
completely fixes from a microscopic level the parameters 
used in the two mode description. 

First let us consider the simplest full GP simulation, 
listed as I in Table H] This consists of the three com- 
ponents starting from the same initial population imbal- 
ances and basically gives a similar Josephson tunneling 
behavior for the three components. As can be seen in 
Fig. [H the Josephson regime is fully identified on the cou- 
pled behavior of Zm and S(f>m- Together with the Joseph- 
son oscillation there is a transfer of population between 
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FIG. 2: Full simulation of Eqs. ^ and two-mode analysis of 
some cases listed in Table |l] The first/second, third/fourth 
and fifth/sixth rows correspond to simulations Ila, Ilia, IVa, 
respectively. Solid lines correspond to the GP simulations. 
Dashed lines depict two-mode results with the parameters 
computed microscopically as described in the text. In most 
cases the two lines in each panel are almost indistinguishable. 



the three different states, see panel (b) of Fig. [TJ As 
discussed above the population transfer dynamics decou- 
ples from the Josephson tunnehng in this regime and thus 
allows to clearly identify the value of NU2, which is of 
course directly linked to C2. The agreement between the 
two-mode and the full GP simulation is remarkable as 
can be seen in Fig. [TJ Taking into account that for *^Rb 
|c2| << Co and therefore U2N « UqN, it is easy to 
prove from the above two-mode equations that, for this 
case, the behavior of the imbalance of all the compo- 
nents follows: Zm — —tOjZm with UJj = UJr\/l + NUq / UJr ■ 

Which corresponds to the Josephson frequency of a scalar 
condensate completely decoupled from the population 
transfer Q. Therefore, the Josephson tunneling is di- 
rectly related to the spin independent coupling, propor- 
tional to Uq. 

Now we consider three distinct cases: Ila, Illa, IVa. 
As listed in Table U they correspond to different initial 
population imbalances for the three components and to 
a different initial number of atoms populating each sub- 
level from the one used in I. In figure we show the 



TABLE I: Conditions of the different full spinor GP simula- 
tions, Eqs. ([!]). (5(jim(0) = in all cases. 



Sim 


iVo(0)/7V 


^-i(0) 


20 (0) 


2+1(0) 


Transfer 


I 


0.4 


0.005 


0.005 


0.005 


YES 


Ila(b) 
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YES(NO) 
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0.000 
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FIG. 3: The first/second, and third/fourth rows correspond 
to simulations Ila(b) and IVa(b) described in Table |T] Solid 
(black) lines correspond to Ila and IVa while dashed (red) 
lines stand for lib and IVb, which do not include the popula- 
tion transfer terms. 

results of the full GP simulations (solid lines). Rmis Ila 
and Illa produce essentially Josephson tunneling dynam- 
ics modulated by a longer oscillation. Simulation IVa, 
produces a much longer tunneling, the ±1 components 
remain mostly on their original side of the trap while the 
one remains mostly balanced. In the first two cases the 
oscillations of the phase differences are fully characterized 
by CO J. In the same figure, and almost indistinguishable 
from the full GP results, we present the predictions of 
the two-mode model. 

As mentioned above the population transfer dynam- 
ics fully decouples from the Josephson tunneling of the 
three components in the considered conditions. Its coun- 
terpart is however not true, the Josephson dynamics gets 
affected by the population transfer as we will discuss in 
the following. 

To clearly see the effect of the population transfer 
terms on top of the Josephson tunneling dynamics we 
consider the same configurations, labeled as "a", but 
without the population transfer terms, "b". The two- 
mode model, without the corresponding transfer terms, 
also reproduces the dynamics of the "b" runs. In Fig. [3] 
we depict in all cases a comparison between the full GP 
solution and the same case but neglecting the population 
transfer term. 

The effects of population transfer are clearly seen on 
the evolution of z„i. In simulation II, which has 2:0(0) = 
it is observed that the long oscillation which modulates 
the full runs, ujti is not present when we switch off the 
transfer term. Instead the population imbalance shows a 
Josephson-like tunneling oscillation which for i ^ 100 ms 
looses the small Zm regime. Therefore, the transfer term 
tends to stabilize the Josephson-like behavior over longer 
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FIG. 4; Figure which shows the frequencies entering in the 
problem, (a) Time evolution of the number of atoms popu- 
lating the m = — 1 sublevel in simulation I, solid line. The 
dashed line depicts a cos(ti;Ti) which is the two-mode predic- 
tion for A'o ~ N/2, (b) Full GP evolution for z_i of simula- 
tion I, solid line. The dashed line shows a cos(ajjt) behavior, 
clearly identifying the Josephson time scale, (c) Full GP evo- 
lution of z-i for simulation II. The dynamics is governed by 
{ujT,iL}j)- (d) The solid black (red-dashed) line corresponds 
to the GP evolution of z-i (z+i) of simulation IV. The dotted 
lines follow a cos{u!rt), which drives the long-time scale of the 
problem. The scales in the vertical axes are not shown for 
clarity. 

periods of time. The absence of transfer of populations 
does not show up on the behavior of the phase difference, 
as can be seen in Fig. [3J which mostly follows the same 
evolution as for the GP equations with the transfer term. 

As in the case of binary mixtures , taking opposite 
initial imbalances for the m = ±1 components enhances 



the Rabi like oscillation and cancels the Josephson one. 
Simulation IV corresponds to such a case, with z_i(0) = 
— z+i(0) and zo{0) = 0. The Rabi oscillation gives rise 
to a long tunneling behavior but in this case modulated 
by the lot oscillation, as can be seen in Fig. [3] and in the 
lowest panel of Fig. 21 If we switch off the transfer term 
the lut oscillation disappears and the limit of small z and 
6(1) becomes unstable. 

Finally, Fig. 3] summarizes the relevant frequencies 
which enter in the interplay between Josephson tunnel- 
ing and population transfer dynamics in the considered 
regime. The first panel isolates lot = NU2, governing the 
transfer of populations, whereas the second one shows wj, 
which sets the fast behavior of the imbalances. The third 
panel shows z^i from simulation II, which is dominated 
by {ijJT,i^j) and the fourth one shows both z±i from sim- 
ulation IV, that are dominated by two frequencies {lot, 

LOr). 

The ability to perform an experiment with spinor 
F — 1 BEG in the conditions considered in this work 
would present for the first time the combined effects of 
Josephson tunneling phenomena and the transfer of pop- 
ulation between different Zeeman components of a spinor 
condensate: the decoupling of the exchange of popula- 
tions from the Josephson dynamics, the identification of 
the different time scales and the role of the population 
transfer in the stability of the Josephson oscillations. In 
addition, a precise measurement of the population imbal- 
ances and global populations of the three species would 
provide an alternative access to the microscopic proper- 
ties of the atom-atom interactions. 
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